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ABSTRACT 

In the past decade, surveys of the stellar component of the Galaxy have revealed a number of streams 
from tidally disrupted dwarf galaxies and globular clusters. Simulations of hierarchical structure 
formation in ACDM cosmologies predict that the dark matter halo of a galaxy like the Milky Way 
contains hundreds of subhalos with masses of ~ 10® Mq and greater, and it has been suggested that 
the existence of coherent tidal streams is incompatible with the expected abundance of substructure. 
We investigate the effects of dark matter substructure on tidal streams by simulating the disruption 
of a self-gravitating satellite on a wide range of orbits in different host models both with and without 
substructure. We find that the halo shape and the specific orbital path more strongly determine the 
overall degree of disruption of the satellite than does the presence or absence of substructure, i.e., the 
changes in the large-scale properties of the tidal debris due to substructure are small compared to 
variations in the debris from different orbits in a smooth potential. Substructure typically leads to an 
increase in the degree of dumpiness of the tidal debris in sky projection, and in some cases a more 
compact distribution in line-of-sight velocity. Substructure also leads to differences in the location 
of sections of debris compared to the results of the smooth halo model, which may have important 
implications for the interpretation of observed tidal streams. A unique signature of the presence of 
substructure in the halo which may be detectable by upcoming surveys is identified. We conclude, 
however, that predicted levels of substructure are consistent with a detection of a coherent tidal stream 
from a dwarf galaxy. 

Subject headings: cosmology: theory — dark matter — galaxies: structure — methods: N-body simu- 
lations 

1. INTRODUCTION 

In recent years, surveys of the stellar component of the 
Milky Way have reveale d rich structure and a number 
of new tidal streams ( Bclo kurov eFa T 20061 1 iGrillmaii] 
2006t [Grillmair fc Dion atosI l2006b : Belokurov et al. 
2007a| ). In addition to the well-studied Sagittarius 
CSgr) stream (llbata et all 12 001': 'Maiewski et al." '2003|; 
iNewberg et all 120031: iMar tfncz-Dclgado et al. 2001 
and the Monoce r os stream (jNew bcrg ct al. 2002|; 
lYanny et al.l l2003t iPenarrubia et al.. 2005.) . streams 
from the tidal disruption of globular clusters such 
as P al omar 5 (|Odenkirchen et al] [2001': 'Rockosi et aLl 
2002t lOdenkirchen et all l2003t iGr illmair & Dionatoi 
2006al) and NGC 54 66 (jGrillmair fc .Tohnsoni 12001 
Belokurov et al.|[2006aD have also been identified. Tidal 
streams probe the gravitational potential on large scales, 
and consequently may prove to be a valuable tool for 
constraining the mass distribution of the Galaxy. 

Simulations of structure formation in ACDM cosmolo- 
gies predict as many as several hundred dense clumps 
with masses > 10® M© in the halo, roughly an order of 
magnitude more than the observed abundance of dwarf 
galaxies surround ing the Milky Way (jKlvpin et al.lll999l : 
iMoore et al.lll999h . One class of solutions to this 'missing 
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satellites' problem involves modifications of the proper- 
ties of the dark matter partic le, e.g., warm dark mat- 
ter (iHogan fc Dalcantonl I200C ) or self-interacting dark 



matter ( Spergel fc Steinhardtl 



2000), in order to reduce 



the amount of substructure. Suppression in the small- 
scale end of the primordial power spectrum has also 
been i nvoked to reduce the predicted abundance of sub- 
halos (jKamionkowski fc Liddlil2000l : IZentner fc Bullockl 
l2003f ). Alternatively, mechanisms by which only a 
fraction of the dark matter subhalos have visible stel- 
lar components at the present e poch have been pro- 
posed to explain thi s discrepancy (jBullock et al.l[200lU 
iKravtsov et"aLll2004D . In these scenarios, a large number 
of dark substructures are predicted to populate the halo. 

In the past few years the number of known satellites 
of the Milky Way has increased dramatically with the 
discovery of 11 new, faint objects whic h are likely to 
be dw ar f galaxies or globula r clusters (IWillman et al.l 
^OSaira; IZucker et al.l l 2006al |bl: iBelokurov et all l2006d . 
2007b; Irwin ct al. 2007), hinting that a significant num- 
ber of very faint or dark satellites could populate the 
halo. If such a population does exist, these structures 
may be detectable in several ways. Gamma rays from 
the annihilation of dark matter particles i n substructures 
may reveal their presence^ ( e.g., Bcrgstrc Tm et al.l [19991: 
'Caicaneo-Roldan & Moore' 2000; Tasitsiomi fc Olintd 
2002; Stochr c t al.^ 2003; Koushiappas ct al. 200i 
iDiemand et al I l2007t l. Dark substructures may also 
be detectable via flux an omalies in mult i- image 
gravitational lens systems (iMao fc Schneider Il998t 



iMetc alf fc Madaul 120011: iDalal fc Kochanek"2002^. 'Chibal 
2002.') or via spectroscopic lensing (Moustakas fc MctcaQ 



2003( 1. However, recent work favors microlensing by stars 
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in the lens galaxy over dark m atter substructure as an 
explanation for flux anomalies (jPoolev et al.ll2067l ). 

In this study we take a theoretical approach to deter- 
mine whether substructure as predicted by ACDM mod- 
els may be detectable by its influence on the formation 
and properties of tidal streams from dwarf galaxies. We 
select a variety of orbits in several host galaxy models 
and simulate the tidal disruption of a progenitor satellite 
on these orbits both with and without dark matter sub- 
structure to evaluate the influences of the orbital path, 
host potential, and presence of substructure on the tidal 
debris. In f|2]we discuss the phase space properties of or- 
bits in the context of tidal disruption, describe our host 
galaxy models, and outline our method for selecting a 
diverse range of orbits. We give the details of our nu- 
merical simulations in and present our results in 21 
We discuss our study in relation to previous work on this 
subject in [J5l and finally summarize our conclusions in 

m 



2.1. 



2. ORBITS 
Phase space properties of orbits 



In this study we use some basic principles governing 
the motion of orbits in phase space to select a range of 
orbits which may lead to streams with diverse properties. 
The structure of phase space in Hamiltonian dynamics is 
well studied, and we refer the reader to some of the many 
excel l ent texts on this subject (e.g.. iBinnev fc Tremaind 
119871: iLichtenberg fc LiebermanHlOQa T for a more thor- 
ough review than is presented here. 

In a spherically symmetric potential, orbits conserve 
four integrals of motion (the energy, as well as each of 
the three components of angular momentum), so orbits 
in such potentials are generically rosettes confined to a 
single plane in configuration space. In a non-spherical 
potential most orbits conserve three integrals of mo- 
tion {regular orbits) or fewer {chaotic orbits), and thus 
densely fill 3-dimensional volumes in configuration space. 

A resonant orbit is a regular orbit whose fundamen- 
tal oscillation frequencies in configuration space are re- 
lated via a linear equation such as Iwi + muj2 + nuj^ = 0, 
where l,m,n are integers. Such a relation between the 
fundamental frequencies of motion implies that only two 
of the three frequencies are independent, and is simi- 
lar to requiring that a fourth integral of motion be con- 
served. Consequently, a resonant orbit satisfying one or 
two resonance equations as above will be confined to a 
sheet or a closed lo op, re s pectively, in configurati on space 
(iValluri fc Merrittl[l998HMerritt fc Vallurilil999D . 

In most potentials that are near-integrable, the vast 
majority of orbits are regular orbits. Only a tiny frac- 
tion of phase space is expecte d to be occupied by chaotic 
orbits and resonant orbits (|Lichtenberg fc Liebermanl 
fl99l) . More numerous than the resonant or chaotic 
orbits, however, are the regular orbits which are 'reso- 
nantly trapped'. These orbits are confined to an island 
in phase space around a resonance and typically occupy 
only thin volumes in configuration space, and so tidal 
debris from satellites on these orbits would be similarly 
confined. Furthermore, if a satellite's orbit is associated 
with a stable resonance, the tidal debris along the orbit 
may be less susceptible to perturbations by substructure 
than streams formed along non-resonant orbits. Theoret- 
ical work has shown that the density of streams on reso- 



nant o rbits decreases more slow ly than on non-resonant 
orbits (jVogelsberger et al.l[2008h . 

Despite the differences in the phase space proper- 
ties of chaotic and regular orbits, it is generally not 
possible to distinguish between these classes without 
full orbit information for >, 25 dynamical time s (e.g., 
iPapaphilippou fc Laska^ 119961: IValluri fc Merrittj flQQS). 
far longer than the expected age of detected tidal 
streams. For this reason, we do not distinguish be- 
tween chaotic and regular orbits, and instead consider 
only whether a given orbit is resonantly trapped. 

2.2. Orbit selection 

We constr uct libraries of orbit s using a method simi- 
lar to that of' Cretton et al.l ([1999) , and select from these 
libraries orbits on which to initialize the progenitor satel- 
lite in our numerical simulations. Each library consists 
of orbits with a fixed energy E and ^-component of an- 
gular momentum L^. Fixing these two quantities in a 
given potential determines a set of coordinates which are 
turning points for orbits (where w_r = = 0), termed 
the zero velocity curve (ZVC). A library is generated by 
initializing orbits from various positions on the ZVC, and 
integrating these orbits in the host potential. 

Our host potential consists of a smooth NFW 

dark matter halo (Navarro et all Il997f ). a Miyamoto- 
Nagai disk (iMivamoto fc Nagailll975f ). and a Hernquist 
bulge ()Hernauistlll990f) . The masses of these components 
and halo concentration are adopted from the favored 
model of iKlypin et al.l ()2002i ) which gives the halo mass 
A/^haio = 10""^^ Mq and concentration c = 12, the disk 
mass Md = 4 x 10^ ° M0, and the bulge mass {mi+m2 in 
iKlvpin et alll2002[ ) Mf, = 8 x 10^ M©. The scale lengths 
of the Miyamoto-Nagai disk, a = 6.5 kpc and b = 0.26 
kpc, and the Hernq uist bulge sc ale param eter, a = 0.7 
kpc, a re taken from i Johnston et a l. (1996i). iKlvpin et al.l 
()2002D use exponential density profiles for the baryonic 
components of the Galaxy, but since these models are 
considerably more computationally expensive to imple- 
ment in our simulations, we instead use the Miyamoto- 
Nagai disk and Hernquist bulge to describe the baryon 
distribution. These parameters produce a rotation curve 
w hich is qualitatively similar, but not identical, to that 
of lKlvpin et all (12002 ). 

Tidal debris can be dispersed both by precession of 
the orbit of the progenitor due to a non-spherical potcn- 
tial and by heating from interactions with substructure 
(jibata et al.ll2002l : iMaver et al.ll2002f ). so it is necessary 
to evaluate the relative importance of these two influ- 
ences. ACDM simulations of structure formation flnd 
that galaxy-sized dark matter halos are generally triax- 



^^^^^^^^^^ 

pendent (e.g.. 


Frenk et al.ll988l:lKazantzidis et al.ll2004al: 


iHavashi et al. 


I2OO7D. We consider the simpler case of an 



axisymmetric dark matter halo with the axis ratio inde- 
pendent of radius. The lengths of the principal axes are 
deflned to he {x : y : z) — {a : b : c) with the x and 
y axes equal and the ^-coordinate aligned with the axis 
of rotation of the Galaxy. The axis ratio q = c/a = c/b 
is varied to generate orbit libraries for halos with q — \ 
(spherical), q = 1.3 (prolate), and q = 0.6 (oblate). The 
potentials for the oblate and prolate versions of the NFW 
halo are obtained by deflning a variable f which replaces 
r in the spherically-symmetric NFW potential 'I'nfw {"f) i 
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TABLE 1 
Orbit Parameters 



Halo Orbit rapo »'pcri SOS 
q (kpc) (kpc) 




0.6 


Ol 


86.2 


23.2 


NR 




02 


74.7 


32.8 


R 




03 


77.8 


30.1 


NR 




04 


54.2 


43.2 


R 


1.0 


SI 


96.6 


15.2 


NR 




S2 


74.2 


44.9 


NR 




S3 


61.4 


58.5 


R 


1.3 


PI 


99.0 


15.2 


NR 




P2 


102.1 


18.7 


NR 




P3 


101.7 


20.6 


R 




P4 


102.9 


23.2 


R 



where P = + y'^ + {zjqf . 

The energy of an orbit is set by equating it to the 
energy of a circular orbit in the plane of the disk at a 
radius of Rem- In this work we focus on orbits with the 
energy of a circular orbit with i?circ = 60 kpc and = 
0.5Lj, max, which produces a range of orbits with diverse 
properties, and integrate 100 orbits in each library. 

From these libraries we select orbits on which to ini- 
tialize our progenitor usin g surface of section plots (e.g., 
iBinnev fc Tremaind Il987l ). Since resonant orbits con- 
serve an additional integral of motion, their path col- 
lapses to a single point (or a small number of points, for 
higher-order resonances) in a surface of section. Reso- 
nantly trapped orbits occupy a small region surrounding 
the parent orbit in the surface of section. The surface of 
section for the orbits in our oblate halo model library is 
shown in Figure [TJ The two concentric ellipses centered 
at Wfl.sos = in the surface of section are generated by 
orbits that are resonantly trapped with an orbit at the 
center of these ellipses. This halo model also admits an- 
other resonant orbit family at roughly i?sos — 45 kpc, 
Wfl.sos = 100 km/s. 

Properties of the orbits we have selected to simulate 
are given in Table [T] with orbits which appear to be res- 
onantly trapped labeled R in the surface of section (SOS) 
column, and non-resonant orbits labeled NR. Also given 
are the apocenter (rapo) and pericenter (r-pori) distances 
of each orbit. Projections of the orbital path and sur- 
face of section for some individual orbits are shown in 
Figure H 

While our method does not exactly identify any or 
all resonances in a given potential, we expect that 
the orbits we identify as resonantly trapped are suffi- 
ciently close to phase space resonances that they will ex- 
hibit the resonant behavior discussed in ^ 32. II Previous 
work has shown that relatively sma ll perturbations do 
not d estroy the strongest resonances (jSideris fc Kandrupl 
|2004[ ). therefore even though we choose orbital initial 
conditions from surfaces of section in smooth static po- 
tentials, we are confident that when we add substructure 
the resonantly trapped orbits continue to be associated 
with the same resonances, although their proximity to 
the resonance may change. 

3. NUMERICAL SIMULATIONS 

We simulate the tidal disruption of our satellites usin g 
the N-body simulation code GADGET-2 (|Springe]||2005[ ) . 



Fig. 1. — Surface of section for 100 orbits in a model with a 
dark matter halo with q = 0.6. All orbits have been integrated for 
at least 50 orbital periods. A large family of resonantly trapped 
orbits is seen in a 'resonant island', which appears as concentric 
ellipses centered on {Rbob,vr^sos) — (45, 100). 

which calculates gravitational forces using a hierarchical 
multipole expansion algorithm. The host galaxy model 
described in ^2.2\ is incorporated with a static potential 
using a modification to GADGET-2 kindly provided by 
B. Robertson. 

The initial satellite is modeled with 5 x 10^ parti- 
cles following a NFW density profile with virial mass 
Afvir — 10^° Mq and concentration parameter c = 12. 
Our procedure for bringing the initial satellite into quasi- 
equilibrium with the host potential, described below, re- 
sults in a progenitor of ~ 3 x 10^ Mq, roughly the mass 
of a dwarf galaxy such as Sgr. 

In this work we focus on a dwarf galaxy as our pro- 
genitor, as oppos ed to a globular cluster as in, e.g., 
Ilbata et all (|2002f ). because we expect that a dwarf 
galaxy will produce debris which is more spatially ex- 
tended due to its substantially greater mass and intrin- 
sic velocity dispersion, and therefore will present a much 
larger cross section for interactions with substructures. 
Since the most massive subhalos are rare and our sim- 
ulations are evolved for only a few orbital periods, this 
increased cross section for interaction may easily be the 
difference between a massive substructure encountering 
the debris at least once or not at all. 

The mass of the NFW profile diverges as r oo, 
so we truncate the satellite's density profile outside of 
the virial radius r yir usin g an exponential cutoff as in 
iKazantzidis et ah! (|2004b( ) . We set the radial scale which 
determines the sharpness of the cutoff rdccay = 0.35rvir, 
which results in a total mass of ^ 1.1 x Afvir- The 
satellite initial conditions are generated using t he ex- 
act di stribution function as detailed in Kazantzidis et aTl 
(|2004bf l. and we acknowledge S. Kazantzidis for provid- 
ing us with code which implements this method. The 
gravitational softening length for the satellite particles is 
set to Esat = 0.07 kpc. 

This initialization scheme sets up a satellite which is 
in equilibrium when isolated, but will not be in equilib- 
rium when placed in the host potential. As a result, the 
satellite will initially lose mass at an artificially high rate 
due to the sudden introduction of an external tidal field 
until it has relaxed into the host potential. To minimize 
this effect, we create a progenitor in quasi-equilibrium 
with the host potential. We place the initial satellite on 
circular orbits in the plane of the disk at 30, 40, and 50 
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Fig. 2. — Projections in x-y, x-z, r-z, and surface of section i?sos-i'ij,sos for selected orbits from each halo model. SI is a regular rosette 
orbit and appears as a nicely filled curve which covers a large range of radii and velocities in the SOS plot. Orbit S3 is close to a resonance 
and occupies a much smaller region in the SOS plot (note that the axes have been rescaled). 02 is extremely close to the strong resonance 
seen in the upper part of Figure [T] and 03 is an orbit that lies at the edge of the resonant island. Orbit P3 is close to a resonance and 
appears as three points in the SOS plot, while non-resonant orbit PI fills a large area of the SOS plot. 



kpc, and track the mass loss rate over time. For all three 
cases there is an initial period of rapid mass loss, but the 
mass loss rate stabilizes after less than 0.2 orbital peri- 
ods. We generate a 'remnant' to use as our progenitor by 
selecting all particles still bound to the satellite after 0.2 
orbital periods for the satellite on the 50 kpc orbit. The 
closest distance at which we initialize any of our orbits 
is ^ 54 kpc, so we choose the 50 kpc orbit for generat- 
ing the remnant to temper the most intense initialization 
shock our progenitor would experience. The remant we 



use as our progenitor contains ^ 140k particles and has 
a mass of ~ 3 x 10^ Mq. 

For our simulations in halos with substructure, wc in- 
corporate as softened point masses 253 subhalos from 
halo G2 in the cos mological N-body simulation of 
iKravtsov et al.l (|2004[ ) and thank A. Kravtsov for mak- 
ing this data available to us. The mass range of these 
subhalos is ~ lO^'-lO^o M©. Since our host potential is 
not identical to that of the host in the cosmological sim- 
ulation (e.g., ours includes baryonic components), using 
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TABLE 2 

SuBHALO Group Properties 



Mass Range Number of e^uh lisoftlmax |aNPw|max 

(IO^OMo) Subhalos (kpc) (10^ km/s/Gyr) {10^ km/s/Gyr) 



4.62 1 6.64 3.03 3.07 

0.1 - 1.16 12 3.57 2.61 2.64 

0.01 - 0.1 67 1.25 1.85 1.88 

0.0008 - 0.01 173 0.46 1.38 1.40 



the final positions and velocities of the subhalos in the 
cosmological simulation as our initial conditions will not 
produce identical orbits. Instead, we maintain the final 
distribution of the subhalos in configuration space and in 
the z-component of angular momentum L^: we initialize 
each subhalo in our simulations at its final spatial coor- 
dinates in the cosmological simulation, and transform its 
velocity by requiring that the current position is a turn- 
ing point of the orbit (vr = Vz — 0) and that its Lz is 

p reserved. 

iNavarro et al.l (|1997D showed that density profiles of 
CDM halos are well-f it by a NFW profi l e over a wide 
range of masses and iKazantzidis et all (|2004bf ) found 
that dark matter subhalos retain their steep inner density 
profile even when tidally stripped. Thus ideally we would 
model the subhalos in our simulation using a NFW den- 
sity profile, but for numerical simplicity we model the 
subhalos as spline-softened point masses. The subhalo 
particles are sorted by mass to form four groups with 
different softenings. The subhalo softening length esuh 
for each group is chosen such that the maximum accel- 
eration due to the most massive spline-softened subhalo 
losoftlmax in each group is equal to the NFW accelera- 
tion at r = O-Olr^ (a radius well-wit hin that shown to 
be sta ble against tidal disruption by IKazantzidis et all 
I20n4hl ^. 

The spline-softened acceleration is maximized at r ~ 
0.8e, while the NFW acceleration increases toward a fi- 
nite limit as r —^ 0, so using a spline softened poten- 
tial increases the effective size of the subhalos. We note 
that for the NFW profile the acceleration at r = O.Olrs 
is ~ 99% of the maximum achieved by the NFW pro- 
file loNFwImax- Howcvcr, large accelerations are reached 
even at moderate radii for a NFW profile: e.g., the NFW 
acceleration attains ~ 40% of the maximum at r = Tj. 
The functional form of the softened potential requires the 
softened acceleration to exceed the NFW acceleration for 
a small range of radii in order for the softened acceler- 
ation to anywhere reach the value for the NFW model 
at r = O.Olrg. Requiring the softened acceleration to be 
smaller than the NFW acceleration at all radii would re- 
strict the maximum softened acceleration to significantly 
less than the NFW maximum, so in order to more fully 
investigate the effect that the steep inner profile of NFW 
subhalos may have on tidal streams we allow the softened 
acceleration to exceed the NFW value at some radii. 

The softening for each group is set by first determining 
the concentration parameter c for the most massive sub- 
ha lo in the group using the c(Mvir) relation determined 
bv lBullock eraIT(|2001ah . which gives 18 < c < 40 for our 
subhalos. The NFW scale radius is then determined 
from c and M^ir for the most massive subhalo, and the 
softening length is then set from as described above. 



Table [2] summarizes the characteristics of each group. 

We do not explicitly include a stellar component in 
our satellite model, but we attempt to follow the distri- 
bution of a stellar component by identifying a subset of 
dark matter particles as 'star particles'. In the initial 
equilibrium halo, star particles are tagged according to a 
Hernquist distribution. After the remnant is generated 
there are ~12,000 star particles remaining. 

For each selected orbit, we simulate the disruption of 
the remnant progenitor in the host model both with and 
without dark matter subhalos, evolving our simulations 
for 4.9 Gyr. 

4. RESULTS 

4.1. Sky distribution of the debris 

The projection of the star particles in Galactic coor- 
dinates after 4.9 Gyr of evolution for a subset of orbits 
is shown in Figures [3] (spherical halo model) , |4] (oblate 
halo model), and [5] (prolate halo model). Comparing the 
simulations in all three halo models, it is immediately 
notable that even those without subhalos do not con- 
sistently produce thin, coherent streams. Simulations 
with subhalos often lead to clumpier structures in these 
coordinates, but the degree of additional scattering or 
clumping for the same orbit simulated with and without 
subhalos is small compared to the variation seen between 
different orbits and different halo models without subha- 
los. 

In Galactic coordinates, the S3 orbit without subhalos 
produces the most coherent stream, which is not surpris- 
ing given that this orbit has rapo ''-^ ?'peri 60 kpc. The 
SI orbit is largely confined to the plane of the disk of 
the Galaxy and has a large apocenter to pericenter ratio 
(?'apo/?'pori ^ 6), and in this projection produces clumpy 
debris clustered along the Galactic plane in the smooth 
halo model. For both orbits in the spherical potential the 
presence of substructure produces some additional dis- 
persion. Substructure also leads to noticeably clumpier 
debris for the S3 orbit, and results in debris displaced 
from the Galactic plane for orbit SI. 

In the oblate and prolate halo models, the resulting 
streams are quite dispersed compared to the S3 stream 
and somewhat clumpy both in the simulations with and 
without subhalos. For the orbits shown in the oblate 
halo, 02 and 03, as well as the PI orbit in the pro- 
late halo, the debris appears more clustered in the sim- 
ulations with subhalos. Orbit P3 produces somewhat 
smooth streams which are qualitatively very similar both 
with and without subhalos, suggesting it has not had any 
significant interactions with the subhalo population. 

Resonantly trapped orbits S3 and P3 produce slightly 
more coherent streams than the non-resonant orbits in 
the same halos. Yet, both orbit 02, which is closely 
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Fig. 3. — Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in the spherical halo model. Orbits 
in models with subhalos often produce clumpier debris than in models with a smooth halo component. The additional disruption due to 
substructure is of a similar degree for the resonantly trapped (S3) and non-resonant (SI) orbits. The SI orbit with subhalos has clearly 
been perturbed, displacing a large amount of the debris relative to the smooth halo simulation. 





Fig. 4. — Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in the oblate halo model. The debris 
for orbits in models with subhalos appears slightly more clustered in this projection than for the models without subhalos, with most 
of the debris following a single path. Both the resonajitly trapped orbit (02) and the non-resonant orbit (03) produce debris that is 
fairly dispersed, regardless of the presence of subhalos. The 03 simulation with subhalos produces an interesting feature at ^ ~ —45°, 
0° < 6 < 30° resembling a bifurcation in the stream. 



7 



+90 



-ISO 




-180 ;, 




+90 




P3 



+1S0 



+90 



180 : 




P3 SUBS 



+180 



Fig. 5. — Sky projection in Galactic coordinates of star particles after 4.9 Gyr for selected orbits in the prolate halo model. The addition 
of subhalos results in slightly more clustered debris in the case of PI, while the debris from orbit P3 qualitatively shows little difference 
between models with and without subhalos, producing streams in both cases. 



associated with a resonance, and the non-resonant orbit 
03 produce quahtatively similar debris which is fairly 
dispersed. In this limited sample of orbits we do not 
find that resonantly trapped orbits consistently produce 
substantially more coherent streams than non-resonant 
orbits, and it does not appear that these orbits are less 
susceptible to the effects of substructure than the non- 
resonant orbits. 

In many cases, orbits in models with subhalos result 
in sections of debris in very different locations than in 
the smooth halo models. The SI orbit with subhalos 
produces a stream that is clearly offset from the Galac- 
tic plane for £ > 0°, while the stream from the same 
orbit in the model without substructure wraps around 
the Galaxy at 5 ~ 0°. Similar formations can be seen 
in most of the simulations with substructure, suggest- 
ing that disturbances in tidal streams due to interactions 
with substructures (both dark and luminous) in the halo 
may be an important consideration when constructing 
models to fit observations of tidal streams. Of particular 
interest is a small stream of stars slightly offset from the 
main stream at £ 45°, 0° < b < 30° in the 03 sim- 
ulation with subhalos which in these coordinates quali- 
tatively resembles the bifurcation seen in the Sgr stream 
([Fellhauer et al.l [200 61. This feature is not apparent in 
the smooth halo simulation. We do not in any way sug- 
gest that our orbit 03 is an appropriate description of the 
orbit of Sagittarius, but draw attention to this feature as 
an example of the type of effect substructure can have on 
tidal streams which may complicate the interpretation of 
observed tidal streams. 

4.2. Structure in phase space coordinates 

The distribution of the star particles in line-of-sight 
velocity in the rest frame of the Galaxy vios rest and he- 
liocentric distance rgun is shown in Figures 1 [3 and El 
It is notable that simulations with subhalos tend to pro- 
duce debris which is less coherent in these coordinates 
than debris in their corresponding simulations in the 
smooth halo models. However, the variation in the pat- 



terns formed in the absence of substructure is quite large 
compared to the changes induced by the addition of sub- 
structure, as was the case for the sky distribution of the 
debris. For example, orbits SI and S3 in the smooth 
halo models produce drastically different patterns: the 
SI simulation produces clear bands of debris in these co- 
ordinates, while the S3 simulation appears as two com- 
pact clumps with little internal structure. 

Many simulations with subhalos are distinguishable in 
these coordinates by a conspicuous trail of debris at large 
radii strongly correlated in uios.rost- Furthermore, uios,rcst 
for the particles in this feature is positive and increases 
with increasing radius, unlike patterns formed by the 
wrapping of debris in smooth halo models (e.g., in the SI 
and P3 orbits without subhalos). This feature is present 
to some degree in most of our simulations with subhalos 
in all halo shapes, and is consistently absent when the 
satellite is evolved in a smooth halo. These 'arms' of star 
particles typically extend far beyond the debris seen in 
the smooth halo simulations. Many orbits in simulations 
with substructure produce debris extending to ~ 200 
kpc, while debris from these orbits in smooth halo simu- 
lations is typically confined to ~ 100 kpc. The particles 
in these arms are generally not very spatially clustered, 
so their strong correlation in velocity may prove essential 
for identifying stars as originating from the same progen- 
itor, as well as for cleanly identifying this feature as the 
result of substructure. We emphasize that this observ- 
able signature appears to be unique to simulations with 
subhalos, making a detection of such correlated stars a 
potential smoking gun indicating the presence of sub- 
structure in the halo, regardless of the shape of the dark 
matter halo. 

It is important to note that some of the dark matter 
subhalos in our simulation may host a luminous compo- 
nent, making them analogous to known dwarf galaxies or 
other satellites of the Milky Way, so this signature could 
be evidence of interactions wit h luminous substruc tures 
as well as dark ones. However. IStrigari et al.l ()2007| ) pro- 
pose that the well-known luminous satellites of the Milky 
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Fig. 6. — Heliocentric line-of-sight velocity in the rest frame of the Galaxy versus heliocentric distance of star particles in the spherical 
halo model after 4.9 Gyr for the same orbits as in Figure [3] Simulations with subhalos produce arms of debris in these coordinates at 
relatively large radii (rsun > 150 kpc) strongly correlated in velocity. This feature is also seen for many orbits in oblate and prolate 
models with subhalos. The S3 simulation with subhalos produces multiple arms of debris, the longest of which extends to ~ 200 kpc. For 
comparison, the debris from the same orbit in the simulation without subhalos scarcely extends beyond ~ 100 kpc. 





Fig. 7. — Heliocentric line-of-sight velocity in the rest frame of the Galaxy versus heliocentric distance of star particles in the oblate halo 
model after 4.9 Gyr for the same orbits as in Figure |4] As in the other halo models, arms of star particles at large radii are present for 
orbits in simulations with subhalos, with the 02 arm stretching past 250 kpc. Neither of these orbits shows significant structure in these 
coordinates when simulated in a smooth halo, even though 02 is a resonantly trapped orbit. 



Way do not all inhabit the most massive dark subha- 
los, which in turn implies that some dark subhalos have 
masses comparable to the masses of the known luminous 
satellites. Under this assumption, any signatures of sub- 
structure detected would not necessarily be more likely 
due to luminous structures than to dark ones. 

To confirm that this striking feature does not result 
from lack of precision of the integration scheme, we re- 
peated the S3 simulation with smaller timesteps to pre- 
vent unphysical two-body interactions, and found the 
same feature present. Additionally, we repeated a num- 



be r of simulations using the subhalos taken from host G3 
in iKravtsov et al.l ()2004f ) , and again found that many of 
the simulations developed similar trails of stars at large 
radii, suggesting that the feature does not arise from a 
rare, serendipitous arrangement of the subhalos. 

To identify the mechanism leading to this feature, we 
examine the acceleration produced by the subhalos on 
selected star particles. We estimate the maximum ac- 
celeration produced by a subhalo on a given star par- 
ticle over the course of the simulation by determining 
the closest approach of the particle to the subhalo in any 
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Fig. 8. — Heliocentric linc-of-sight velocity in the rest frame of the Galaxy versus heliocentric distance of star particles in the prolate 
halo model after 4.9 Gyr for the same orbits as in Figure[5] Like the simulations in oblate and spherical halo models with subhalos, the PI 
simulation with subhalos produces debris at large radii correlated in velocity, but this feature is not seen as clearly in the P3 simulation 
with subhalos. 

arms had a strong encounter with a very massive subhalo 
(> 10^0 M0), while the particle not located in an arm 
(marked with a blue triangle) did not have any encoun- 
ters producing accelerations greater than 300 km/s 
/Gyr recorded in a snapshot. 

The time evolution of the total energy of selected par- 
ticles is shown in Figure I10| along with the proximity of 
each particle to the subhalo which caused it the greatest 
acceleration. The particles selected correspond to those 
in Figure [9l For the three particles in the arms at the 
end of the simulation (the two top panels and the bottom 
left panel), a large increase in energy is evident, clearly 
corresponding to the time of closest approach of the sub- 
halo. The energy of the particle in the bottom right 
panel, which had a lesser maximum acceleration in Fig- 
ure[Hl shows smaller deviations, but which still appear to 
be correlated with the proximity of the subhalo. The net 
change in total energy for this particle over the course of 
the simulation is small compared to the increases in the 
other particles' energies. For the three particles which 
did gain a substantial amount of energy, the majority of 
the increase occurred during a single interaction with the 
subhalo producing the largest recorded acceleration. 

Figure [TT] shows the distribution of the fractional 
change in total energy over the course of the simulation 
for all satellite particles (both dark matter and star par- 
ticles) in the 03 simulation with subhalos. The solid line 
shows the change in energy for all particles, and the red 
dashed line represents only the particles in the longer 
arm at the end of the simulation {vsun > 170 kpc). The 
small positive offset from zero seen for the distribution 
for all particles shows that satellite particles on average 
gained a small amount of energy {'^ 8%), indicating that 
interactions with the subhalo population have a slight 
tendency to increase the energy of the stream particles. 
The right panel is rescaled to show the distribution of 
the particles in the arm in relation to all the satellite 
particles. The arm particles appear to naturally fill the 



Fig. 9. — {color online) Maximum acceleration due to each sub- 
halo experienced by four selected star particles versus subhalo mass 
for the 03 simulation. Particles shown are selected from features 
seen in Figure [7] the shorter, upper arm (black crosses), the tip 
of the long arm {red diamonds), the middle of the long arm {pur- 
ple squares), and the main region {blue triangles). Final rsun and 
iiltjs rest I plotted in Figure [7] are labeled for each particle. 

snapshot. Due to the limited time resolution of our snap- 
shots, the distance of closest approach determined in this 
way is likely to be significantly overestimated, and hence 
the maximum acceleration may be significantly underes- 
timated. 

Figure [5] shows the maximum acceleration due to each 
subhalo as a function of subhalo mass for four selected 
star particles in the 03 simulation with substructure. 
Each of the four particles was selected from a different 
region of the rsun-i'ios,rcst plot in Figure [71 a particle 
from the shorter, upper arm (black crosses), the tip of 
the longer arm (red diamonds), the middle of the longer 
arm (purple squares), and the lower part of the main re- 
gion (blue triangles). A clear trend is evident, with the 
most massive subhalos typically producing the largest 
accelerations. A large acceleration does not guarantee a 
substantial energy exchange, but indicates the strength 
of an individual encounter. All three particles from the 
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Fig. 10. — Total energy per unit mass (solid line) and distance to subhalo producing the largest acceleration (dashed line) for four 
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Fig. 11. — Fractional change in total particle energy for all particles in the 03 simulation with substructure over the course of the 
simulation (solid line) and for all particles in the longer arm at the end of the simulation (red dashed line, right panel, color online). 
Particles defined to be in the arm have rnun > 170 kpc, and include dark matter particles and star particles. Eq is each particle's total 
energy at the beginning of the simulation. The particles in the arm occupy the upper end of this distribution, accounting for the majority 
of particles with an increase in energy of at least ~ 45%. 

positive end of the AE/\Eo\ distribution, accounting for 
almost all of the particles which experienced an energy 
increase of ^ 45% or more. In addition, all of the 55 star 
particles in this feature experienced an acceleration from 
a single subhalo of at least 600 km/s /Gyr, and for all 
but one of these particles the most massive subhalo was 
responsible for the strongest recorded encounter. 

A number of mechanisms for transferring large 
amounts of energy to stream particles are possible, in- 
cluding two-body encounters (i.e., the acceleration and 
energy gain of the particle is due primarily to the influ- 
ence of a single subhalo), three-body encounters (i.e., the 
energy gain is due to two subhalos producing comparable 
accelerations on the particle during the time the major- 
ity of the energy transfer occurs), and less easily cat- 
egorized interactions (e.g., a particle-subhalo encounter 



during which the acceleration of the subhalo by the host 
potential provides a non-negligible contribution to the 
particle's energy gain). We have examined in detail the 
circumstances of the large energy gain for a few selected 
particles, deferring a thorough investigation of the en- 
ergy exchange mechanism to a future study. In our sim- 
ulations we find examples of an isolated two-body inter- 
action between a particle and a subhalo and of a particle- 
subhalo interaction during which the acceleration of the 
host potential appears to play an important role. We do 
not find clear evidence of three-body interactions, but 
our limited inquiry does not rule out the possibility of 
such interactions. 

To test the sensitivity of our results to the subhalo soft- 
ening scheme, we consider the change in a particle's total 
velocity by an isolated interaction with a massive subhalo 
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Fig. 12. — Change in velocity |A-!/| = If/ — for a particle 
scattered by the NFW (solid) and softened (dashed) potentials for 
the most massive subhalo in our simulations. The thick, medium, 
and thin lines are for b = 5, 15, and 25 kpc respectively. 

modeled by a NFW or a softened potential. The angular 
deflection of a particle's trajectory by a central force can 
be determined by integrating the orbit equation: 

e = 7r-2 , (1) 

where E and £ are the total energy and angular momen- 
tum of the particle per unit mass, U(r) is the potential 
per unit mass, and rmin is the distance of closest approach 
of the particle. Defining the impact parameter b to be 
the perpendicular distance between the undeflected path 
of the particle and the center of force, i = h\plE — bvrci 
with Wroi the relative velocity before the encounter. The 
scattering of a particle by a substantially more massive 
subhalo is well-described by this setup as long as the 
subhalo does not accelerate significantly during the in- 
teraction. 

The magnitude of the change in a particle's velocity 
during this interaction \Av\ = ju/ — w^j is determined for 
a given b and v^ci ■ Figure [T^] shows | Av \ as a function 
of Wioi for selected values of b using the NFW and the 
softened potentials for the most massive subhalo in our 
simulations. For an isolated scattering interaction with 
this subhalo, the maximum |Aw| is ~ 60 km/s for the 
NFW potential (at Vrd 110 km/s, & 22 kpc), with 
a fairly weak dependence on Vj-d and b. For the softened 
potential, |At/| is more strongly dependent on Vi-d and 
b, achieving a maximum of ~ 145 km/s (at v^ci 130 
km/s, & ~ 14 kpc), with a typical | Aw | somewhat smaller 
depending on the parameters. \Av\ reaches a minimum 
of 35 km/s for the NFW potential and ~ 43 km/s for 
the softened potential for 5 < 5 < 20 kpc and 100 < 
Vrei < 300 km/s. 

This calculation suggests that our choice of softening 
may lead to changes in the particle's velocity larger by 
up to a factor of ~ 2.5 for the strongest isolated inter- 
actions than if we had implemented the NFW potential. 
However, we emphasize that this calculation only gives 
an estimate of the change in a particle's velocity, which 
cannot be easily translated to an effect on the properties 
of the debris without detailed information about the pa- 
rameters of a particular interaction in the rest frame of 
the Galaxy. The change in the particle's kinetic energy in 
the rest frame of the Galaxy depends on the initial veloc- 
ity vector of the particle in the Galactic rest frame, and 



so these scattering events can reduce as well as increase 
the total energy of the particles. 

From this analytical estimate, we expect that any fea- 
tures in our simulations attributable to interactions with 
subhalos, such as the arms of debris in Figures[6l[7l andUl 
could be somewhat less prominent for subhalos modeled 
with a NFW potential. Because the parameters deter- 
mining the energy gain or loss of a particle may vary 
greatly between scenarios in which scattering events oc- 
cur, it is not possible to determine, in general, how much 
less prominent these features will be for NFW subhalos, 
or how often the resulting debris will differ substantially 
from that obtained using softened subhalos. 

Interactions in which the acceleration of the subhalo 
by the host galaxy contributes appreciably to the energy 
gain of the particle are also seen in our simulations. In 
these cases the two-body scattering approximation is not 
an adequate description of the interaction, and our use 
of a spline-softened instead of a NFW potential for the 
subhalo plays a smaller role in determining the energy 
gain of the particle. 

Previous w ork (Ibata et al.ll2002l; I Johnston et al.ll2002t 
iMaver et al.ll2002: Penarrubia et al.ll2006f l has suggested 
that interactions with substructure would increase the 
velocity dispersion of the streams. To determine whether 
the presence of substructure can be inferred from this 
measurement, we examine the spread in wios.rest as a func- 
tion of Galactic longitude for star particles after 4.9 Gyr, 
shown in Figure [13] for three selected orbits. The S3 or- 
bit, which has 7'apo/''peri 1 and is close to a phase space 
resonance, produces a stream with a velocity spread of 
roughly 100 km/s in the absence of substructure. When 
substructure is added, the stellar debris is projected as 
two clumps with a slightly larger range of velocities in 
these coordinates. Debris from the 02 orbit is, sur- 
prisingly, more compact in velocity when subhalos are 
present than in the smooth halo. Similarly, the P3 orbit 
develops streams which are slightly narrower in velocity 
in the simulation with substructure than in the simu- 
lation without substructure. The PI orbit (not shown) 
produces stellar debris dispersed over ~ 400 km/s in ve- 
locity at all longitudes regardless of the presence of sub- 
structure. Thus a consistent trend towards larger veloc- 
ity dispersions when substructure is present is not seen. 

5. DISCUSSION 

The question of whether tidal streams can be used to 
constrain substructure in the Milky Way has been ad- 
dressed by a number of studies for a variety of scenarios. 
We discuss here the results of some of those works in 
relation to our own. 

The effects of halo shape and substructure on the tidal 
disrup tion of a globular cluster are studied bv llbata et al.l 
(|2002r ). who find that heating due to a population of 
subhalos decouples from the effects of precession in a 
non-spherical potential by uniquely increasing dispersion 
in the z-component of angular momentum. Our study, 
while more sensitive to the effects of isolated interactions 
with massive substruct ures, is not as sen sitive to this 
effect since, as noted by llbata et al.l ()2002[ ). debris from 
a dwarf galaxy is already fairly dispersed in Lz and thus 
the effect of substructures on this quantity is more subtle. 
Those authors conclude that the detection of a single 
cold stream from a globular cluster would provide strong 
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Fig. 13. — Line-of-sight velocity in the rest frame of the Galaxy versus Galactic longitude of star particles after 4.9 Gyr for S3, 02, 
and P3 simulations with and without substructure. Contrary to expectations, narrower streams in these coordinates are produced in the 
presence of substructure than in the smooth halo for orbits 02 and P3. 



evidence against the substructure picture. In contrast, 
by examining a large range of scenarios we find that in 
some cases tidal streams from a dwarf galaxy in a halo 
with substructure remain as cold as streams in a smooth 
halo. However, due to the differences in our approach 
and choice of progenitor satellite, it is unclear whether 
this result is generalizable to streams from a globular 
cluster, and thus further study is necessary to clarify the 
si gnificance of detection s of cold tidal streams. 

iJohnston et al.l (|2002l ) examine the effects of substruc- 
ture on test particles in a live host halo, and find that 
they can distinguish between smooth halos and those 
with substructure using their proposed 'scattering index'. 
Their study examines streams arising from a spherical 
host potential and thus does not compare the dispersion 
caused by substructure with that due to precession, so 
additional tests are necessary to determine the robust- 
ness of the scattering index for uniquely constraining 
substructure. Of particular importance for this study, 
however, is their conclusion that neglecting the response 
of the halo to the substructures, as we have done, is likely 
to underestimate the influence of the substructures due 
to the wake generated by their motion in the halo. 

The properties of tid al streams in an evolving host po- 
tential are examined bv lPeharrubia et all (|2006D . includ- 
ing a study of the effect of substructure on streams in a 
host potential with a spherical dark matter halo. They 
conclude that substructure contributes only modestly to 



heating in , and only see this effect in the trailing arms 
of streams unbound more than 5 Gyr ago, which were 
subjected to a higher rate of encounters with subhalos 
than more recently unbound material. Our implementa- 
tion of substructure does not account for the evolution of 
the subhalo population, and so the frequency of encoun- 
ters does not have an externally set time-dependence, 
but depends rather on the particular orbital path. It 
is demonstrated in Figure [9] that over 5 Gyr a parti- 
cle in our simulation will have interacted to some degree 
with a large number of the subhalos, and consequently 
we assume that our star particles will have collectively 
experienced all but exceptionally rare types of interac- 
tions. Thus we expect that integrating the streams for 
longer times will simply enhance the features due to sub- 
structure that are already apparent. 

6. CONCLUSIONS 

We have investigated the effect of halo substructure 
on the tidal debris from a dwarf galaxy. Our conclusions 
can be summarized as follows: 

1. Our results suggest that the shape of the host po- 
tential and the path of the orbit in question have 
a stronger influence on the overall spatial distri- 
bution of tidal streams than does the presence of 
substructure in the halo. In light of the enormous 
amount of variation in the properties of the tidal 
debris produced by simulations in a smooth halo, it 
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is clear that dispersal of the debris is not an effect 
that can be attributed exclusively, or even primar- 
ily, to substructure. We find that the detection of a 
coherent tidal stream from a dwarf galaxy does not 
rule out the abundance of substructure considered 
in this work. 

2. Substructure can noticeably shift the location of 
sections of debris relative to the debris in the cor- 
responding smooth halo simulation. This is likely 
to be an important complication for constructing 
models to fit observational data. 

3. Interactions with a population of subhalos can lead 
to significant changes in the small-scale properties 
of the tidal debris. Satellites on some orbits de- 
velop rich structure in the Tsun-^^ios.rest distribution 
in a smooth halo which is largely smeared out when 
simulated with substructure. The presence of sub- 
structure also tends to lead to more clustered de- 
bris in a sky projection, and in some cases smaller 
velocity dispersions. 

4. Stars at large distances with a small range of posi- 
tive line-of-sight velocities with respect to the rest 

frame of the Galaxy which increase with increasing 
distance may be a unique signature of an interac- 
tion with a relatively massive halo substructure. 
We find that particles in this configuration typi- 
cally have experienced a single strong encounter 
with a subhalo of mass > 10^° Mq. The formation 
of this feature thus depends on the orbits of the 
most massive subhalos in relation to the spatial and 
velocity distribution of stellar material. The large 
number of tidal streams expected to be present in 
a halo like that of the Milky Way, coupled with the 



fact that stars from the stellar halo of the Galaxy 
would similarly be susceptible to such interactions, 
implies a high likelihood for producing these fea- 
tures. Despite the fact that these particles are ob- 
served at large radial velocities and distances, they 
are not unbound and will thus eventually return. 
Each individual feature is therefore transient. 

5. While satellites on resonantly trapped orbits ap- 
pear to produce more coherent tidal debris than 
non-resonant orbits in some individual cases, we 
do not consistently see this trend in our simulated 
debris. Similarly, we do not find that streams from 
satellites on resonantly trapped orbits are less sus- 
ceptible to the effects of substructure. 

With the wealth of data expected from upcoming mis- 
sions, we are optimistic about the prospects of using tidal 
streams to gain a better understanding of the mass dis- 
tribution of the Milky Way. If the halo of the Galaxy 
contains a population of substructure as predicted by 
ACDM, clear signatures of its presence may be detectable 
in the near future. 
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